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■ We discuss the current situation concerning measurement and readout of Josephson-junction based 
qubits. In particular we focus attention of dispersive low-dissipation techniques involving reflection 

, of radiation from an oscillator circuit coupled to a qubit, allowing single-shot determination of the 

■ state of the qubit. In particular we develop a formalism describing a charge qubit read out by 
Pj^ ' measuring its effective (quantum) capacitance. To exemplify, we also give explicit formulas for the 

readout time. 

(N 

I. INTRODUCTION 

Nanotechnology is considered promising for fabrication of scalable solid-state electronics for quantum computersiS^. 
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However, progress towards solid-state quantum computing will critically depend on the development of measurement 
schemes and readout devices that, on demand, can determine the state of individual qubits in a fraction of the 
coherence time, but which otherwise do not disturb the qubit system. In quantum optics, efficient measurement 
techniques have been developed during the last thirty years based on laser-atom interactions and recently implemented 
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in e.g. ion traps^Si&S. 



The corresponding work for solid state systems effectively started only about ten years ago, and is currently 
exploring various paths. A problem is that there is no general device for operation and readout, like a laser, but 
rather a multitude of implementations of measurements of charge, spin, magnetic flux and charge current that must 
be adapted to the specific qubits to be studied. Therefore, the qubit readout technology must be developed in 
intimate connection with the qubits for characterization and control of coherence properties. This is a painstakingly 
slow process which, however, cannot be circumvented, because it is essential in many respects. In particular, it is an 
l— "~ ', important tool for determining the coherence properties of the qubits. Moreover, the technology not only concerns 
I ' qubit readout devices, but also involves quantum oscillators for storing and transmitting information and for coupling 
, qubits. 

I/"} • Interestingly enough, quantum-optical methods are now being applied to solid-state qubit systems, using microwaves 
OO \ for operating and reading out qubits, and oscillator circuits and transmission lines for coupling qubits, introducing 
cavity-QED in solid-state systemsi2iii*i£. This may turn out to be a major road on the Road Map for quantum 
coherent systems ("quantum computers" and will be at the focus of the present paper. In particular we will 

describe some practical schemes for reflecting microwaves from an oscillator circuit, the phase shift measuring the 
changes in charg o 15 i 16 or magnetic flusiLiSiiSiSSiSi induced by a qubit, allowing to distinguish between the different 
states |0) and |1) of the qubit. 



II. MEASUREMENT OF QUANTUM INFORMATION AND QUBIT READOUT 

A. Introduction 



The ultimate objective of a qubit readout device is to distinguish the eigenstates of a qubit in a single measurement 
"without destroying the qubit", a so called "single-shot" quantum non-demolition (QND) projective measurement. 
This objective is essential for several reasons: state preparation for computation, readout for error correction during 
the calculation, and readout of results at the end of the calculation. Strictly speaking, the QND property is only 
needed if the qubit must be left in an eigenstate after the readout. In a broader sense, readout of a specific qubit 
must of course not demolish other qubits in the system. 

Note that one cannot "read out the state of a qubit" in a single measurement - this is prohibited by quantum 
mechanics. The effect of a single ideal measurement on a qubit in a general superposition a |0) + b |1) is to leave the 
qubit in one of the states |0) or |1), which carries no information about the amplitudes a or b. It takes repeated 
measurements on a large number of replicas of the quantum state to characterize the state of the qubit - "quantum 
tomography" . This is the procedure to collect the statistics for expectation values. 

The measurement connects the qubit with the open system of the detector, which collapses the combined system of 
qubit and measurement device to one of its common eigenstates. If the coupling between the qubit and the detector 
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is weak, the eigenstates are approximately those of the qubit. In general however, one must consider the eigenstates 
of the total qubit-detector system and manipulate gate voltages and fluxes such that the readout measurement is 
performed in a convenient energy eigenbasis (see e.g. Refs. I22I and I23I) . 



B. Survey of readout methods for JJ-based qubits 

Here we will provide a brief recapitulation of the "history" of readout of Joscphson-junction (J J) based quantum 
circuits and qubits. For an extensive discussion of JJ-based qubit circuits, see the recent review by Wendin and 
Shumeiko^. 

Figure 1 shows general designs for the charge and flux qubits and with oscillator-type readout circuits. The Single 
Cooper-pair Box (SCB) (similar to Fig. 1, left) is described by the Hamiltonian 

H = E c {n~n g ) 2 - Ejcos<f> (1) 
and the rf-SQUID (similar to the flux qubit in Fig. 1, right; see Section VI) by the Hamiltonian 

H = E c n 2 — Ej cos + E L ^ - ; (2) 

where Ec is the charging energy of the SCB island, Ej is the Josephson energy due to Cooper-pair tunneling between 
the superconducting electrode and island across the JJ, and El is the inductive energy of the superconducting loop, n 
and <f> are operators for the induced charge and the phase of the (effective) Josephson junction, and n g is the induced 
charge controlled by the external bias voltage V g . Both qubits are represented by the generic 2-level Hamiltonian 

H = --(ea z +Aa x ) (3) 

where a z and a x are the usual Pauli matrices. 




FIG. 1: Circuit diagrams and 2-level energy spectrum of two basic JJ-qubit designs: the SCB charge qubit with LC-oscillator 
readout (left), and persistent-current flux qubit with SQUID oscillator readout (right). For the charge qubit, the control 
variable e on the horizontal axis of the energy spectrum (middle) represents the external gate voltage (induced charge), and 
the splitting is given by the Josephson tunneling energy mixing the charge states. For the flux qubit, the variable e represents 
the external magnetic flux. In both cases, the energy of the qubit can be "tuned" and the working point controlled. Away 
from the origin (asymptotically) the levels represent pure charge states (zero |0) or one |1) Cooper pair on the SCB island) or 
pure flux states (left |0) or right |1) rotating currents in the SQUID ring). 

In the original experiment of Nakamura et al^, demonstrating coherent oscillation of the charge qubit 2-level system 
(Fig. 1, left), the readout was implemented simply by a control dc-pulse on the charge gate, moving the working point 
far away from the origin so that the upper |1) level ended up above the gap edge of a superconducting lead connected 
to the SCB island via a tunnel junction. As a result, a Cooper pair on the upper level |1) would immediately decay 
into the external lead as two quasiparticles, creating a normal electron current. Repeating the measurement at a 
high rate created a detectable current proportional to the occupation of the upper state |1), revealing the oscillations. 
Since the SCB is permanently connected to the environment via a tunnel junction, it seemed plausible at the time 
that this might be the reason for the short coherence time, ~ 2-3 ns. 

This focussed the interest on more advanced readout devices. A remedy could be to use a charge measuring 
device that was only capacitively coupled to the SCB island and could be turned on and off by an external voltage 
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pulse. Delsing and coworkers^ 5 , therefore developed an rf-SET (radio frequency single-electron transistor) readout*^ 
for the charge qubit, and successfully detected free oscillations and studied the detailed behaviour of relaxation and 
dephasing 2 ^. The result showed that the coherence time was confined to below 10 ns and seemed limited by relaxation 
effects. Moreover, subsequent experiments by the NEC group^l, implementing more advanced readout concepts, 
storing the emitted pair of quasiparticles on a superconducting island, and reading the charge with a superconducting 
SET, made no big change. All in all, the status seems to be that the coherence time of the circuit is severely limited 
by intrinsic charge fluctuation processes (noise) in the substrate, or in the tunnel barriers, or by transients due to the 
pulsed operation of the qubit. 

Alternatively one could create a new type of charge qubit by connecting the Cooper Pair Box island to two JJ tunnel 
junctions, creating a Single Cooper Pair Transistor (SCT). This could be probed via charge 16 or curren1i22i22i22i2i 
measurements. These experiments can be designed either as threshold detection measurements or as microwave- 
reflection measurements. The reflection measurements with phase-shift detection will be the main theme of this 
paper. 

Moreover, there is the persistent-current flux qubit, based on a quantum version of the RF-SQUIP 33 34-35 36 coupled 
to a measurement dc-SQUID. This measurement SQUID be operated either as a current threshold detectos 3 ^ 6 - or 
as resonance circuit reflecting and phase-shifting microwave radiation 2 ^!. We will briefly describe the microwave 
reflection measurement also in this case in sect. VI. 



III. CHARGE MEASUREMENTS 



The most straightforward way to read out a charge qubit is to measure its charge. As discussed above, to obtain a 
high fidelity read out one should perform a measurement in the qubit eigenbasis. This removes the possibility for the 
qubit to switch its state during the measurement. When the measurement basis is fixed, as determined by a charge 
measurement, we need to bias the qubit where the charge basis is the eigenbasis. For the superconducting charge qubit 
this implies a complete quenching of the Josephson energy, while for quantum dot charge qubits one needs raise the 
tunnel barrier between the dots. Having quenched the transitions between the charge states the fidelity will in theory 
be perfect, irrespectively of the measurement speed. In reality there is always some remaining transition/relaxation 
channel open which implies the need of a fast read-out. Fast read-out is also mandatory for implementing an error 
correcting algorithm, where the read out and correction should be performed on a time-scale set by decoherence of 
the other qubits. 



A. The radio-frequency single-electron transistor 

The state of the art charge measurement device is the radio-frequency single-electron transistor (rf-SET)^ 6 . with an 
experimentally measured sensitivity of Sq = 3.2 • 1CT 6 e/VlIz 3 ^. The measurement time t ms needed to separate two 
states with a charge difference AQ is t ms = (25q/ AQ) 2 , indicating the possibility of detecting a charge difference of 
one percent of the electron charge (AQ = 0.01 e) in half a microsecond. 




FIG. 2: A single-electron transistor: a small metallic island connected to source and drain leads through tunnel junctions. The 
signal voltage V s induces a charge q m on the measurement capacitance C m . When the system is at the limit of being Coulomb 
blockaded, a small change of q m will have a large effect on the current 7 through the SET. 

The single-electron transistor (SET) consists of a small metallic island connected to source and drain leads through 
tunnel junctions (see Fig. 2). Applying a source-drain voltage the current through the SET depends critically on the 
charge induced on the gate capacitance. The charge on the gate capacitance is determined by measuring the current. 
The SET is in itself a sensitive electrometer but suffers from low operating speed, which in addition to being a 
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FIG. 3: Resonant circuits for read-out: a) A lumped element LC-oscillator coupled to a driving source and a radio-frequency 
detector through a transmission line, b) The radio-frequency single-electron transistor measuring the charge of a charge qubit 
(SCB). The current through the SET determines the dissipation in the resonant circuit. The dissipation is determined by 
measuring the amplitude of the reflected signal, c) Setup for measuring the quantum capacitance of the charge qubit. The 
qubit capacitance influences the resonance frequency of the oscillator. The capacitance is measured by determining the phase- 
shift of the reflected signal. 

drawback on its own makes it sensitive to the low frequency charge fluctuations. The rf-SET is realized by embedding 
an SET in a resonant circuit (see Fig. Fig. 3 a,b). The main source of dissipation in the oscillator is current flowing 
through the SET, and for small amplitude oscillations we may replace the SET by its effective (differential) resistance. 
The oscillator is excited by sending down a radio-frequency signal on resonance. The dissipation is determined by 
measuring the amplitude of the reflected signal. This way of operating the SET increases its operating speed and 
sensitivity significantly. 

B. Single-shot read-out 

The high sensitivity can be used to couple the rf-SET weakly to the charge qubit, reducing the back-action in the 
off-state. In practice it is impossible to switch off the interaction completely, i.e. the qubit eigenbasis is not exactly 
the charge basis, and there is some unavoidable mixing of the charge states . A careful investigation of the rf-SET 
coupled to a superconducting charge qubit shows that single-shot read-out with a very high fidelity is still possible in 
practice^2i22il£. To optimize the fidelity one should bias the charge qubit as far away from the degeneracy point as 
possible, in order to minimize the effect of any residual Josephson coupling. 
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IV. CAPACITANCE MEASUREMENTS - A QUASICLASSICAL DESCRIPTION 

The low-frequency charge fluctuations present in all realizations of superconducting charge qubits so fai4i^£ pro- 
motes the use of schemes where the charge qubits are operated at the charge degeneracy point around the origin in 
Fig. 1 (middle panel). Here the qubit eigenstates have equal average charge and thus they are shielded from charge 
fluctuations. To use rf-SET read-out described in the previous section one needs to quickly shift the qubit far away 
from degeneracy simultaneously with switching-on the measurement. Although this timing is far from impossible, 
the present trend is to use schemes where the qubit remains at charge degeneracy also during read-out in order to 
minimize decoherence. This can be achieved by measuring the effective capacitance of the charge qubit, as we describe 
below. 



A. Quantum capacitance of a single-Cooper-pair box 

The quantum capacitance of the Cooper-pair bo:s^iii is related to the parametric capacitance of small Josephson 
junctions^^S^i which is a dual to the Josephson inductance. The origin of the quantum capacitance of a single- 
Cooper-pair box (SCB) can be understood as follows. Assume that we put a constant voltage V m on the measurement 
capacitance of the SCB, i.e. we put a voltage source between the open circles in Fig. [3]:. The amount of charge on 
the measurement capacitance qm e (V m , V g ) will be a nonlinear function of the voltage V m as well as the gate voltage 
V g and whether the qubit is in the ground or excited state. We may define an effective (differential) capacitance 

C^ e (V m , V g ) = -±- q <>/e(V m ,V g ), (4) 

as seen from the measurement circuitry. Away from the charge degeneracy points of the SCB no charge will float across 
the Josephson junction and the effective capacitance is simply the geometric capacitance C ge0 m = CjC m /{Cj + C m ) of 
the Josephson junction capacitance Cj and the measurement capacitance C m in scries. Around the charge degeneracy 
point a change of voltage will induce a shift of a Cooper-pair across the Josephson junction. For a voltage changing 
slowly on the timescale of the inverse qubit gap fiEJ 1 this charge redistribution is dissipationless. This contribution 
to the effective capacitance, which depends on the qubit state, we call the quantum capacitance Cq S . From the SCB 
Hamiltonian (see e.g. Eq. (|f 8(1 ~) it is straightforward to show 

C^ e (V m , V g ) = ± = Cgoom + C Q (n g + n m ), (5) 

Cj + Cm CS (El[l-2(n g + n m )Y + Ef ' 

where Cs = Cj + C m + C g is the total island capacitance, Eq — 2e 2 /Cs is the Coulomb energy of a Cooper-pair, 
and n g / m = C g / m V g / m /2e are the induced number of Cooper-pairs on the gate and measurements capacitances 
respectively. We note that the quantum capacitance is positive in the ground state and negative in the excited state, 
as illustrated in Fig. 4. The absolute value is largest at the charge degeneracy n g + n m = 0.5 

Cq ~c^W (6) 

The quantum capacitance can be an order of magnitude larger than the geometric capacitance for realistic parameters. 
By inserting the Cooper-pair box in a resonant LC-circuit and detect its influence on the resonance frequency, the 
quantum capacitance was recently measured by two different groups^i^. The measurement setup is similar to the 
setup that would be used for qubit read-out and is analyzed below. 



B. Read-out by measuring the quantum capacitance 

At the charge degeneracy point the effective capacitance of the SCB in the ground and excited state differs by 
2Cg ax . Imbedding the SCB in a resonant circuit as shown in Fig. a) and c) we can detect the corresponding 

change in the oscillators resonance frequency = 1/ \l 1 L(C ± Cg ax ) = =f Cg ax /2C), where cjo = 1/VLC is 

the bare resonance frequency. The voltage reflection amplitude T(lo) — V ou t (w)/Vd(uj) seen from the driving side of 
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FIG. 4: The quantum capacitance of a single-Cooper-pair box with Ej/Eq = 0.2. Here normalized by the maximum value 

Cq -c^eJ- 

the transmission line can for a high quality oscillator be written 

1 + 2iQ^—^- 

r( w ) = . ^° . = e** , where Vr = 2 arctan ( 2Q( " ~ , (7) 

i~2^ ( "'" o) \ ^ J 

up to a constant phase depending on the length of the transmission line. Here Q is the resonator's quality factor, 
which for the circuitry in Fig. 0a) is determined by the characteristic impedance on the transmission line Zq through 
Q = ojqLC 2 /C^Zq. Since there is no dissipation in the oscillator we have |r(w)| = 1. Driving the oscillator at the 
bare resonance frequency uj^ = loq the phase-difference between the ground and excited state of the qubit will be 

8<p r = tp.3 - tp* = 4 arctan (QC^/C) . (8) 

The phase-difference can be detected by measuring the reflected signals in-phase and quadrature components after 
mixing it with the drive. 



V. A QUANTUM DESCRIPTION OF THE QUANTUM CAPACITANCE READOUT 

Above we described the quantum capacitance of the Cooper-pair box and its use for qubit read-out in a quasiclassical 
manner, treating the oscillator, transmission line, drive and detection classically. In order to address questions about 
the optimal read-out time, quantum back-action on the qubit and the quantum efficiency of the read-out process we 
need a fully quantum description of the system. In this paper we do not have the space to go into details, which will 
be published elsewhere^, but we will discuss the principles of our model and show a fully quantum derivation of the 
quantum capacitance. 

The approach we chose is close to the "Quantum Network Theory" introduced by Yurke and Denker— . In sec- 
tion lV Al we start by writing down the Lagrangian describing the classical dynamics of the circuit. Through a Legendre 
transform we arrive at the corresponding Hamiltonian. By stating canonical commutation relations between our phase 
coordinates $i and their canonical conjugate momenta (charges) qi 

[$i,ft]=*ft (9) 

we arrive at a quantum Hamiltonian description of our circuit, which is discussed in section TVBI In the relevant pa- 
rameter regime we arrive at the expression for the quantum capacitance of the Cooper-pair box. Finally in section lV CI 
we give an expression for the optimal qubit read-out time using homodyne detection. 
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A. Circuit Lagrangian 



The circuit for performing read-out through the quantum capacitance is presented in figure A Josephson charge 
qubit is capacitatively coupled to a harmonic oscillator, which is coupled to a transmission line. Through this line, all 
measurement on the qubit is performed. We model the line as a semi-infinite line of LC-circuits in series. The working 
point of the Josephson junction can be chosen using the bias V g . In writing down the Lagrangian we are free to chose 
any quantities as our coordinates as long as they give a full description of our circuit. Since we are treating a system 
including a Josephson junction, the phases $i(t) = / dt'Vi{t') across the circuit elements are natural coordinates, 
as discussed by Devoret in ref. |^£J. (This is in contrast to the original work by Yurke and Denker where charges are 
chosen as coordinates.) 
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FIG. 5: The circuit used for measurement of the quantum capacitance of the Cooper-pair box. It is similar to the circuit 
shown in Fig. |3Ji)+c), but the transmission line is here modelled as a semi-infinite line of LC-circuits in series. The phases 
= J dt'Vi(t') across the different circuit elements are the coordinates used in the Lagrangian describing the dynamics of 
the system. 



The capacitive energy of the circuit act as the kinetic terms in the Lagrangian 3 

T = - ((7 9 $ 2 + Cj&j + C m & m + C<1> 2 + cMn) + 2 E AxC T ($f) 2 



and the inductive part plays the part of potential energy 



V = ^-E J cos( 2 ^ J 



2L T (Ax) 2 

Applying Kirchoff's voltage law gives us the constraints 

* 9 - *j + V g = 0, $j + $ m - $ = 0, 
+ = 0, $ - $ L = 0, 

which gives the Lagrangian for the system 

. _ C qb &, (C ose + C c )$ 2 C C ($P) 2 $ 2 



(10) 



(11) 



(12) 



2e 

2L^ EjC0S [n* J 



(13) 



The capacitances in the Lag rangian are now Cqb — Cj ~\~ Cg ~\~ Cm> and Cose — C -t- Cm,' 

B. Hamiltonian and quantum capacitance 

From the Lagrangian we easily obtain the Hamiltonian through a Legendre transform. We present the Hamiltonian 
on the form 



H — H q b + H, 



Htl + 



(14) 
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and for simplicity we assume weak coupling C m <C {C osc , C q b] and present only the lowest order (in C m /{C q b, C osc }) 
terms. H q b contains the qubit degrees of freedom including the coupling of the qubit to the rest of the system 

H gb = ^(qj + C g V g f + (q + q p ) (qj + C g V g ) - Ej COS (^<Pj) , (15) 

A^qb ^qb^osc \ / 

while the terms describing the oscillator, transmission line and their interaction are 

H ^ ^ H - $ + 1 yf (g (W 2 , (^ +1 -^) 2 \ „ _ qq P (w) 
Hosc -2CZ~ c + 2L> HTL -2C c + Ai^{-^- + 2L T J' Hint ~ c£? (16) 

where the charge operators q, qj, q p and qf are the conjugate momenta to phase operators <E>, <I>j, $^ an d $f respec- 
tively. For realizing a charge qubit the box charging energy Eq = e 2 /2C q b is much smaller than the Josephson energy 
Ec S> Ej. For relevant parameters we can then limit the qubit charge qj to the two values {0, 2}e, and we get the 
usual expression for the qubit Hamiltonian in the language of the Pauli spin matrices 

„ E e i Ej q + q p q + q p 

H qb = — TT a z 7T a x + 2E c k t-a z + 2E c k ^(1 - n ), (17) 

2 2 e e 

where we introduce the electrostatic energy-difference of the qubit states E e i — 4Ec{l — no), the dimensionless charge 
no = C g V g /e, and the oscillator-qubit coupling coefficient k = C m /C osc <C 1. The last term does not influence 
the systems dynamics and may be absorbed in a small shift of q and q p . Rotating the remaining two terms to the 
eigenbasis of the qubit 



H qb = a z ^\QE 2 c {I -no- n(q + q p )/e) 2 + E 2 /4, (18) 

we arrive at the usual charge qubit Hamiltonian with the charge induced by the oscillator added to the induced gate 
charge. We now concentrate on the case when the oscillator frequency ujo = l/y/LC OS c is much lower than qubit 
frequency Ej/fi. Furthermore we consider the amplitude of the oscillators charge oscillations go such that the induced 
charge oscillations in the qubit are small nqo/e <C Ej/AEc < 1. In this case the qubit will follow the oscillator 
dynamics adiabatically and the rates for transition between the qubit eigenstates are negligible. Furthermore, we may 
Taylor expand the qubit energy around the working point and at the charge degeneracy n = 1 we arrive at the final 
Hamiltonian 

H = -- az+ {^7 c + ^E7 a > ( q + qp) + 2L + 2C c + ^{-^ + 2L T 1 (19) 

v ' i=i v * 

The qubit thus shifts the capacitative energy of the oscillator, which in turn corresponds to adding a small extra 
capacitor to the oscillator Cs = C osc + Cq, where the quantum capacitance Cq is given by 

9p 2 C 2 

Cq = -ettw^ (2°) 



which is identical to the semiclassical formula in Eq. |SJ This in turn will shift the resonance frequency of the oscillator 
with an amount 5u> = — <j z ujoCq /2C osc . 



C. Qubit read-out using homodyne detection 

Taking the continuum limit Ax — > in Eq. I|19|) the solutions to the Hamiltonian for the transmission line correspond 
to fields <fr(x ± vt) propagating to the left and right with velocity v = \j \JCtLt- From the Hamiltonian we derive 
equations of motion from which the right-propagating out-field $ out (t) and the charge on the oscillator q(t) can 
expressed in terms of the left-propagating in-fields <£> m (£), which is determined by the drive. In the parameter regime 
relevant for qubit read-out, where the amplitude of the qubit charge oscillations induced by the LC-circuit is small, 
the qubit will follow the oscillator adiabatically. Furthermore it is adequate to neglect third and higher order terms 
in the Taylor expansion of the qubit energy in Eq. Ijl8(l . Thus we arrive at linear equations of motion which we may 
solve in Fourier representation 

i2C c C^Luj 3 
~ 1 - (C s + C c )Loj 2 - iojC c Zo{\ - LCy,oj 2 ) ' 
q(Lo)+q p (Lo) = xM$;"M, 

*"•*(«) = 4rWM = S(Lo)% n (Lo), (21) 
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where Zq — \J Lt / Ct is the characteristic impedance of the transmission line. Since there is no dissipation in the 
lumped circuit we have |<J> "*(u;)| = |$ m (w)|. In this linear regime the Heisenberg equations of motion are similar to 
the classical ones and for a high quality oscillator we have S(uj) = T(uj) as given in Eq. 0. To discuss the quantum 
statistics of the qubit read-out we need a quantized representation of the fields 



<T>"'(/. = ^("^[C^ t + (C) t e iwt ] 



f. iy poo j 



<nn ^| °|[xMCe- iut + ^) , (C) t ^ > (22) 



where the in-field annihilation and creation operators obey the canonical commutation relations 

Kj, a\ji] = S(u - uj 1 ), and [a u , a u >] = 0. (23) 

What we have achieved is a complete quantum description of the dynamics of the oscillator and qubit in terms 
of the incident field from the transmission line. We also get a full description of the outgoing field, which is what 
will enter the detector. The linearity of the equations of motion arise since we have approximated the qubit with a 
state-dependent capacitance. This is valid as long as the amplitude of the charge oscillations induced on the qubit 
island are small, so that a second order Taylor expansion of the energy in Eq. (|18|l is enough. Furthermore we neglect 
transitions between the qubit states, which can be done for small amplitude oscillations and a low oscillator frequency 
huj <C Ej. 

We have now the formalism needed to address questions about the readout time which we do in the next section. 
We can also calculate the backaction of the measurement process on the qubit. This and the issue of the quantum 
efficiency, i. e. the relation between the qubit dephasing rate and the measurement time will be addressed elsewhere*. 

By measuring the power of the reflected signal mixed with a local oscillator, e.g. the drive itself, and then a local 
oscillator shifted 90 degrees the in-phase and quadrature signal amplitude can be extracted, as was done in Ref. I2M 
For optimized qubit read-out it's advantageous to implement the standard quantum optics scheme of homodyne 
detection^. 

To model a measurement we thus put the in-field in a coherent Glauber stated 

\{a(oj)}) = exp (J doj[a(u)(a^ - |0), (24) 

where a(u>) is the Fourier-transform of our drive signal, and |0) is the continuum vacuum field a w |0) = 0. We model 
our drive source with a narrow (Td -C ojd) Gaussian distribution in frequency 

UJd e -{u-u d )^/2Vl 

a(uj) = a - -= , (25) 

i-d V w 

where oiq is a dimensionless constant. For |i| -C T^ 1 this gives the average electrical field 



V m {x,t) = -aoy/2hZaU)dwa.[ud{t + x/v)], (26) 

giving on average T™ = a^LUd photons per second sent through the transmission line by the drive. 

The annihilation operator for the reflected signal has the amplitude a° u * = S{u>)a™ s» e^^ dt ^ v ^ "-'a™ depending on 
the qubit state. The signal is mixed with a strong local oscillator with amplitude a,Loe~^ Jdt ~ VLO ' and the intensity 
is detected. The result is then integrated for a time T. The intensity is given by the number of photons incident on 
the detector 

N{uj,T) = ( dt tf(t)b(t) (27) 



where the field at the detector is 

b{t) = ra LO {t) + rv{t) + ta out {t), (28) 
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here r <C 1 denotes the small reflection coefficient of the mixer and t is the corresponding transmission coefficient, and 
v(t) is the vaccuum-field of the idle mixer port. (The mixer is a beam-splitter in the quantum optics case.) We now 
assume that \raLo(t)\ S> 1 is large and neglect second order contributions in v(t) and a(t) to equation (|2"7j> . arriving 
at the average photon number at the detector 

{N a/e ) = T (r 2 a 2 LO + 2ira LO < cos (tf /e ~ <p L oj) , (29) 

where a™ is amplitude of the in-field. The first term is a pure local oscillator term and contains no information, 
while the second term is maximized choosing ip^o = Vr 1 2 ' Pr — 5. The two probability distributions of the number of 
detected photons will be separated by two times the variance after the measurement time 



T — 



V m A ■ 2 

"4 sin 



2 



(30) 



where T™ is rate of photons sent into the transmission line by the drive. For not too low quality factor of the oscillator 
we can use the Breit-Wigner approximation for x(uS) leading to the expression for the phase-shift in Eq. JHJ). Thus we 
find for the measurement time 



T™ 4 sin 2 [2arctanx] 16T™ C 

From the condition that the drive should induce only small oscillations of the qubit charge we arrive at the following 
bound for the drive strength 

K n <^ + ^ V ) (32) 



giving a lower bound on measurement time 



167T 



T ms > -^(x + x- 1 ), (33) 



E 



indicating that the measurement time must be larger than h/Ej, which is not very restrictive. Due to the oscillator 
"ring-up" time the measurement time is further limited by T ms > Q/luq = xC /CqUq. For a fixed u>q this indicates 
that the regime x <C 1 is advantageous. Comparing these two inequalities we find the shortest measurement time for 



x opt = — Qo P t = xjc^^Ji and T ms > —, (34) 



implying that a low Q is clearly an interesting regime. For low Q the Breit-Wigner approximation of x(w) breaks 
down, and so does the simple estimates of the measurement, but the formalism developed here is still applicable using 
the full expressions. An optimization including the measurement induced back-action on the qubit will be discussed 
in Ref.El 



D. A comparison with dispersive readout using a non-adiabatic oscillator 

In an experiment at Yale university^ 3 - a charge qubit coupled capacitatively to a microstrip cavity was read out by 
sending microwaves through the cavity. The qubit state influences the resonance frequency and thus the phase-shift 
of the transmitted signal. This phase difference was then detected in a similar fashion as described above. 

The main difference compared to what was discussed above is that the cavity resonance frequency (5.4 GHz) was 
of the same order of magnitude as the qubit frequency (4.3 GHz). In this regime it is appropriate to use the rotating 
wave approximation and the system dynamics is described by the Jaynes-Cummings hamiltonianiS. 

For a comparison we start with the Hamiltonian of a qubit coupled transversely to a harmonic oscillator 

H = — o z + iga x [a' - a) + huj osc \a'a+ - I , (35) 
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applicable for a charge qubit at charge degeneracy, coupled capacitatively to the harmonic oscillator. Performing a 
straightforward second order perturbation expansion in the coupling term we find the renormalizcd spectrum 



E n 1 — TibJosc 
E n i = TllOosc 



Tiw, 



gb 



2 



1 



flu,' 



gb 



TujJ r , 



Tilt). 
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qb " 

n 



hui n 



and 



Tiao, 



qb 



tUi)„ 



hui. 



qb 



hui n 



(36) 



where n indicates the number of photons in the oscillator and t / I the qubit in the ground/ excited state for the 
unperturbed state (g — > 0). The spectrum is formed by two equidistant sets of energy levels, where the effective 
qubit-dependent frequency shift of the oscillator amounts to 



Suj 



1 



Huj. 



qb 



flUJ n 



+ 



1 



TlLO, 



gb 



TlUJn 



(37) 



being negative for the qubit in the ground state. In the regime \uj q b — lu osc \ <C 0J q b only the first term contributes and 
we arrive at the Jaynes-Cummings result Suj = —a z g 2 / (huj q b — huj osc ). In the adiabatic regime ll> osc ui q b we can 
neglect the terms ±hu> osc in the denominators giving the result below Eq. (|20() . Suj = —cr z 2g 2 /huj q i ) . Thus we see that 
the frequency shift given by the Jaynes-Cummings hamiltonian can be described on the same footing as the one given 
by the quantum capacitance. 



VI. FLUX MEASUREMENT 



The 2-level quantum states of the persistent-current flux qubit (Fig. 1) are characterized by different directions of 
the persistent currents circulating in the qubit loop, hence different directions of the induced magnetic flux. The flux 
qubit readout is based on the detection of the induced flux or direct measurement of the persistent currents. The 
latter method is also relevant for charge qubits with loop-shape electrodes (e.g. quantronium2&) where the intensity of 
the induced flux is too small to be detectable while detection of the persistent current is possible. In flux qubits with 
larger Josephson junctions, persistent currents are large, and the measurement of flux is not that difficult 18 - 20-35 ' 55 . 
The experimental measurement setup is sketched in Fig. 1, left picture: the qubit loop is inductively coupled to a 
dc-SQUID connected to a current source. Direction of the persistent current in the qubit loop affects the magnetic 
flux threading the SQUID and thus affects the SQUID critical current as well as its plasma frequency. This allows 
one to make the two types of measurements, by probing the dc and the ac properties of the measurement SQUID. 
In the first case, the critical current is measured by applying a dc current slowly increasing with time, and detecting 
the value of the critical current when the SQUID switches to the resistive branch (threshold detection), repeating the 
measurement to create a histogram of the events. In the second case, an ac current is applied and the phase shift of 
the reflected signal is measured. The latter method is also possible to realize using a linear LC-oscillator instead of 
the dc SQUIDS. 




To quantitatively analyze the circuit (see e.g. Refs. I3l l2lll54r) . let us for simplicity consider the single junction flux 
qubit, the analysis also applies to the experimental three junction qubits. The circuit Hamiltonian consists of the 
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Hamiltonian H q of the qubit loop, the Hamiltonian of the SQUID, Hs, and the Hamiltonian of the inductive coupling, 
H int . The qubit Hamiltonian has the form given in Eq. J5J), here we modify the notations, 

H q = E Cq n 2 q - E, Jq cos <p q + E Lq 4> 2 q , (38) 

where Ec q = (2e) 2 /2C q is the charging energy of the qubit junction, Ej q is the Josephson energy, E^ q = 
(<f>o/27r) 2 (l/2i g ) is the inductive energy of the qubit loop, and <fi q — <f> q — f q is the induced flux, f q is the bias- 
ing flux. The SQUID Hamiltonian is given by equation, 

2 ^ 

Hs = J2 ( E C* n i " E Ja COS <A0 + E Ls4>l + ~ <fe), (39) 

i 

where the induced flux is <p s = <j)\ + 4>2 ~ /s> and I(t) is a (non)stationary current bias. The interaction term has the 
form 

Hint = E M 4> q <ps, (40) 

where the interaction energy is determined by the mutual inductance M, Em = {^o/2-k) 2 {M / L q L s ). 

Now we truncate the Hilbert space of the circuit Hamiltonian to the lowest energy states, which include the two 
almost degenerate (for Ec q <C Ej q and at f q ~ tt) lowest energy states in the potential wells of the qubit potential 
energy, Fig. 6, and the ground state of the SQUID. Then the qubit Hamiltonian takes the form, in the eigenbasis of 
the non-coupled wells 



H q = --(ea z + Aa x ). (41) 



Here e(/ g ) is the energy level difference in the wells proportional to f q — tt , and A is the energy level splitting due to 
the macroscopic tunneling between the wells. The truncated interaction term takes the form, 

Hint = E M (/) Q 4> S <7 Z , (42) 

where </>o is the half distance between the minima of the potential energy. The off-diagonal term in the interaction is 
neglected since it is exponentially small. 

To truncate the SQUID Hamiltonian, we introduce new variables, (f>± = (l/2)(0i ± (^2), and <p s = 2<fi+ — f s , 

H s = \E Ca n 2 _ + 2E Cs n 2 s - 2E Js cos 0_ cos + E Ls cj> 2 s + ^I{t)<t>- ■ (43) 

Now we assume that the inductive energy is sufficiently large, E^ s 3> Ej s , to provide small fluctuation of the induced 
flux, (j> s -C 1. This allows us to expand the cosine term; then keeping the first order term with respect to </> s (non- 
vanishing for / s 7^ 0), and taking into account the interaction term l4*2l we write the </> s -dependent part of the total 
Hamiltonian on the form, 

H(4>„) = 2E Cs n 2 s + (e Js cos0_ sin (J^j + E M M.f q )^ + E Ls <f> 2 s , (44) 

This is the linear oscillator shifted from the origin, the shift being proportional to the induced flux in the qubit loop. 
Making projection on the ground state of this Hamiltonian, we arrive at the non-trivial, qubit-dependent part having 
the form, 

H(<t> s )^-2E Js \sin^ a z cos0_, A = ||^. (45) 

Combining this with the rest of the SQUID Hamiltonian and the truncated qubit Hamiltonian, we finally get: 
1 , . , 1 



H = -7,{™z + Acr *) + ^E Cs n± - 2E Js 



(46) 



For the 3-junction qubit^ the coupling constant A in Eq. (|45|l acquires an additional factor— Ejg/E^g, which results 
from tracing out the plasma mode in the qubit loop. This mode does not form the qubit in the 3-junction circuit (in 
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contrast to the single-junction qubit), but this mode is an auxiliary one connecting the qubit to the outside world, 
and it is eliminated similar to the SQUID variable cf> + in Eqs. I|44|) and (|45|l . 

The Hamiltonian (|4fj|> describes a flux qubit directly coupled via an effective coupling constant A to a non-linear 
Josephson oscillator. The coupling affects the Josephson energy of the oscillator; hence the critical bias current, i.e. 
the the magnitude of the bias dc current at which the oscillator switches to the dissipative regime. Quantitatively, 
these critical current values for the 3-junction qubit are, 

Y/c = 1E Js (cos (f s /2) ± sin (/ s /2)) . (47) 

The advantage of this method is that the measurement circuit can be disconnected during the time period between 
the measurements by switching off the flux through the SQUID, f s = 0, thus enhancing the decoherence time of 
the qubit. The disadvantage of the method is that for slow readout (low Josephson plasma frequency of the SQUID 
compared to the qubit frequency) the switching current depends on the average value of the induced flux, (er 2 ), 
which equals zero at the degeneracy point, e = 0. Thus measurement can only be performed by departing from the 
degeneracy point, which is undesirable due to enhanced decoherence. This difficulty can be solved by probing the 
qubit quantum inductance, which is analogous to the quantum capacitance measurement for charge qubits. 

VII. INDUCTANCE MEASUREMENT 




FIG. 7: Single-contact flux qubit inductively coupled to a linear oscillator. 



Consider first a simpler circuit with a linear LC-oscillator replacing the dc SQUID (Fig. 7). Such a device, rf-SQUID 
inductively coupled to a linear oscillator is a classical device employed for many years for precise measurement of 
magnetic fieldS. The principle of operation is based on the magnetic field dependence of the Josephson inductance, 
which affects the resonance frequency of the oscillator probed by an external rf signal. The classical Hamiltonian for 
this circuit, again assuming for simplicity a single Josephson junction in the qubit loop, has the form, 

H = E Cq n\ ~ E Jq cos0 g + E Lq 4> 2 q + E M ^ q <j> + E Cosc n 2 + E Losc (j) 2 + j-J(t)^. (48) 

Neglecting the junction capacitance energy, and expanding the qubit potential energy near the minimum point 4> q o(f q ), 



E Jq cos0 9 + E Lq 4> 2 q = Ufa) = U{4> qQ ) + -U"$ q0 ){4> q - <M 2 , (49) 



we define an effective inductance of the qubit circuit L? q '' via the relation 



-[/"(<M = (<f>o/2vr) 2 /2i^ / (/ 9 ) = E e //(f q ). (50) 



After having diagonalized the total potential energy, we obtain the shift of the oscillator inductive energy due to 
coupling to the qubit, 

E 2 

^° SC ^ LOSC ~4<W (51) 
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This gives rise to a shift of the oscillator resonance frequency depending on the magnetic flux through the qubit loop, 
which is probed with an external rf signal I(t). 

A similar measurement procedure also applies to the quantum regime. The quantum Hamiltonian for the same 
circuit has the form, taking into account Eqs. I|41|) and (|42[1 . 

H = ~(ea z + Aa x ) + E M <j> <t)a z + E C o SC n 2 + E Losc 4> 2 + ^-I{t)<j>. (52) 
2 2e 

For a slow oscillator and weak coupling, the Hamiltonian can be rotated to a qubit eigenbasis, and expanded with 
respect to the coupling term, 

H = -\ea z - {EM2 <f> 2 a z + E Cosc n 2 + E Losc cf + A/(t)0 (53) 
2 e 2e 

(here the term linear in 4> is omitted since it only produces a non-essential small shift of the oscillator coordinate; 
£ = \/e 2 + A 2 ) . The second term in Eq. ll ~> ; > f) provides the shift of the oscillator inductive energy depending on the 
qubit state, 

E 2 (f> 2 

Elosc — Elosc — — O*. (54) 

£ 

Comparing this quantum result with the classical equation (|51|l . we are able to identify the quantum inductance of 
the qubit, 

2tt \ 2 24>l 



Lq = [*- ) -T a *' (55) 

The quantum inductance is inversely proportional to the qubit level splitting, similar to the quantum capacitance of 
charge qubits, and it approaches its maximum value at the degeneracy point. 

This conclusion also applies to the 3-j unction flux qubil^; the only difference is due to the suppressed coupling 
of the qubit to the outside world discussed in the previous section below Eq. (|46|l . In this case, an additional small 
factor, Ej q /EL q , appears in the the coupling term in Eq. (|52|l . giving rise to the following equation for the quantum 
inductance of the 3-junction qubit, 

Proceeding with the case of the readout dc SQUID, Eq. (|46[) . we find that this case seems to be qualitatively 
different from the LC-oscillator readout: the qubit-meter coupling is non- linear. Obviously, this results from the fact 
that the qubit is not directly coupled to the readout </>_-oscillator, but rather via an intermediate cj) s -oscillator, the 
elimination of which results in the non-linear coupling. 

One way to solve the problem is to displace the 4>- -oscillator by applying a constant current bias, I(t) = Iq + Ii(t). 
Then expanding the potential energy around the minimum point, 0_ = <f>- + 0, where (f>- satisfies the equation, 
2Ej s cos(/ s /2) sin0_ = —(h/2e)Io, we arrive, in the linear approximation, at a Hamiltonian similar to Eq. I|52|l . 

H = -\{ea z + Aa x ) + ]-E Cs n 2 _ + E Js cos0_ cos(/ s /2)0 2 + XE Js 9a z + ^h{t)9, (57) 
2 2 2e 



A = A^£itan(/ s /2). (58) 

Another solution would be to access directly the (f> s -oscillator linearly coupled to the qubit. This can be done, for 
example, by driving a bias flux through the SQUID, f s {t). For small variation of phases in Ea. H43|) . the oscillators 
decouple, and the relevant part of the Hamiltonian, taking into account Eqs. H41JI and (j42(l . approaches a form similar 
to Eq. (O, 

H = -hea z + Aa x ) + +2E Cs n 2 s + E Js <f> 2 s + Em ^ j ^ ^ + 2E Js f s {t)cp s . (59) 

2 &Lq 
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Finally, the <p s -oscillator may be accessed by using an asymmetric SQUID, with different inductances of the left and 
right legs (Fig. 1), Li ^ L 2 , L\ + L 2 = L s . In this case, the coupling of the SQUID to an external current source in 
Eq. ffity becomes modified to 

(60) 

which results in direct coupling of the probing current to the <f> s -oscillator. At zero flux bias, f s = 0, the two 
4>± -oscillators do not interact, and the relevant part of the Hamiltonian reads, 

H = -\{ea z + Aa x ) + 2E Cs n 2 s + {E Ls + E Js )<j> 2 s + Em ^ j ^° ^ s(Jz + L * ~ 1 <f, s . (61) 

Z tbLq Ze L s 

Thus there are several ways to employ the dc SQUID for dispersive measurement of the qubit quantum inductance. 
Note, however, that in the latter case the inductive energy of the SQUID plays a role, and for small SQUID inductance 
(El s ^> Ej s ) the oscillator frequency may become large, violating the adiabatic regime assumed in the derivation. 

VIII. CONCLUDING REMARKS 

In this paper we have outlined some practical schemes for capacitive and inductive readout, detecting the state of a 
qubit by reflecting microwaves from an oscillator circuit, the phase shift measuring the changes in charge or magnetic 
flux induced by a qubit, allowing to distinguish between the different states |0) and |1) of the qubit. In particular 
we focussed attention on the single-Cooper-pair box (SCB), the effective capacitance of which can be defined as the 
derivative of the induced charge with respect to gate voltage. In addition to the geometric capacitance, there is the 
quantum capacitance due to the level dispersion at the anti-crossing caused by the Josephson coupling. We described 
the process of reflection of quantized radiation and derived expressions for the shortest measurement time needed to 
resolve the qubit states, suggesting that a low Q should be advantageous for weak back action fast readout. 
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